Cross-sectional STM of Mn-doped GaAs: theory and experiment 
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We report first-principles calculations of the energetics and simulated scanning tunneling mi- 
croscopy (STM) images for Mn dopants near the GaAs (110) surface, and compare the results 
with cross-sectional STM images. The Mn configurations considered here include substitutionals, 
interstitials, and complexes of substitutionals and interstitials in the first three layers near the sur- 
face. Based on detailed comparisons of the simulated and experimental images, we identify three 
types of Mn configurations imaged at the surface: (1) single Mn substitutionals, (2) pairs of Mn 
substitutionals, and (3) complexes of Mn substitutionals and interstitials. 

PACS numbers: 68.37.Ef,75.50.Pp,71.55.Eq 
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I. INTRODUCTION 

GaAs can be doped with Mn to form a dilute mag- 
netic semiconductor with a Curie temperature as high 
as 160 K j^— ir It is generally accepted that holes created 
by the substitution of Mn for Ga mediate the ferromag- 
netic interaction between Mn dopants in this material. ''■^ 
Naively, each substitutional Mn is expected to produce 
one hole; therefore the nominal hole concentration, p, 
should be equal to the Mn concentration. For typical 
Mn fractions of 5%, this implies a hole concentration 
p — 1.1 X 10^^ cm~^. Measured hole concentrations 
of as- grown material are much less than this, typically 
by a factor ^^3. The source of this compensation re- 
mains somewhat controversial, having been variously at- 
tributed to either excess As in the form of antisites and 
interstitialaSiLSia or Mn interstitialsi^°db^^i" Recent ex- 
periments show that careful annealing near the growth 
temperature can significantly enhance the conductivity 
and hole concentrationi^iii This, as well as recent ion 
channelling experimentsjii suggests that interstitial Mn 
is the more likely source. A complete picture is lacking, 
however; for example, the distribution of Mn interstitials 
in as-grown material is still unclear. 

Cross-sectional scanning tunneling microscopy 
(XSTM) is an effective tool for addressing this issue 
because it can image, with atomic resolution, the 
structural and electronic configuration of impurities and 
defects as they are present in the bulk. There are several 
recent XSTM studies of Mn-doped GaAsjiiiii^ these 
studies used the correlation between transport data, 
level of Mn doping, and observed defect density to infer 
the location (substitutional versus interstitial) of Mn 
dopants and As antisites. The relative abundance of As 
antisites and interstitial Mn deduced by these studies 
lacks a clear consensus: one study finds that antisites 
are not present at all,^^- while another finds that both 
antisites and interstitial Mn are present, and compensate 
substitutional Mnii^ 

To investigate the nature of dopants and defects in 
this material, we performed high resolution XSTM mea- 
surements on a (110) cleavage plane of Mn-doped GaAs 
and used first-principles calculations to interpret the im- 



ages. Specifically, we used density-functional theory to 
simulate STM images for a number of Mn configurations 
near the GaAs (110) surface. Configurations were chosen 
based on their calculated energetic stability, and included 
several metastable configurations that would probably be 
kinetically stabilized. A detailed comparison of the re- 
sulting simulated and experimental images leads us to 
the following conclusions: (1) isolated substitutional Mn 
as well as pairs of substitutional Mn occur with roughly 
comparable frequency in bulk Mn-doped GaAs, and (2) 
interstitial Mn is typically bound in complexes with sub- 
stitutional Mn in several different configurations. 



II. EXPERIMENTAL XSTM IMAGES 

The Mn-doped GaAs sample used in this study was 
grown by molecular-beam epitaxy (MBE) using well es- 
tablished methods and conditions iAJ^ A buffer layer of 
high quality undoped GaAs was first grown on an n- 
type GaAs substrate. This was followed by growth of 
a 100 nm thick GaAs buffer layer grown at 250° C and a 
260 nm thick layer of Mn-doped GaAs at the same tem- 
perature. Growth quality was monitored by reflection 
high-energy electron diffraction. The Mn concentration 
was determined with x-ray diffraction to be 0.6%. 

This sample was then characterized by XSTM mea- 
surements in a multichamber ultrahigh-vacuum facility. 
It order to obtain an atomically abrupt cleave across the 
epilayer, the samples were thinned ex situ to <200 /im. 
After being loaded into the XSTM chamber (base pres- 
sure <10~^'' Torr), the samples were scribed and cleaved 
in situ to expose a (110) surface (perpendicular to the 
MBE growth direction), as previously describediJ^ The 
constant-current (40 pA) images shown here are of filled 
electronic states measured with a sample bias of —2.5 V. 

Figure □ shows the results of fiUed-state XSTM mea- 
surements on this sample. The scan area of this image is 
170 X 170 and displays ~1500 surface atoms. Filled 
state STM images of GaAs (110) reveal the As surface 
sublattice; the surface Ga sublattice evident in empty- 
state images, on the other hand, is not directly revealed 
in this image. Since only the top As atoms contribute. 
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FIG. 1: Constant-current fiUed-state XSTM image of a (110) 
cleavage plane of Mn-doped GaAs (001). Five reproducible 
features (A-E) are marked on the image and drawn schemat- 
ically below. 



STM of the GaAs (110) surface reveals only every other 
layer in the (001) direction. 

Figure ^ reveals many surface features not normally 
seen in STM of the GaAs (110) surface. Several occur 
with enough regularity to warrant close attention. Fea- 
ture A is an elliptical region of intensity in an As (001) 
plane with major and minor axes the size of two and one 
As surface atoms, respectively. Feature B is circular, and 
only slightly larger and more intense than the contribu- 
tions from the surrounding As. Feature C is more diffuse 
and broad than both A and B, with intensity spread over 
two pairs of As atoms in neighboring (001) planes. As 
discussed in more detail below, this is what one would 
expect if there were two A features in neighboring As 
(001) planes. Feature D is of similar shape to that of C 
but is more asymmetric, with one side of the defect the 
size of a single surface As atom. There does not appear 
to be a preferred crystalline orientation of feature D as 
different orientations are visible in Fig.^ Feature E is a 
dark region close to a surface As atom displaced in the 
(001) direction; there is also an overall apparent depres- 
sion in the vicinity of this feature, characteristic of band 
bending that occurs near positively charged defects on a 
p-type surfaceiiS 

On the basis of the experimental data alone it is ex- 
tremely difficult, if not impossible, to determine what 
dopants or defects produce these features. As we show in 
this work, however, when interpreted with first-principles 



theory, these measurements resolve both the type and 
crystalline orientation of Mn dopants in GaAs. 

III. THEORETICAL APPROACH 
A. Physical and Electronic Structure 

We modeled the GaAs (110) surface in a supercell slab 
geometry consisting of 5 layers of GaAs in a 4x4 surface 
unit cell. A vacuum region of 13.4 A was used for all 
the calculations. In such a cell the separation between 
a dopant and its periodic image is ~16 A, sufficiently 
large to ensure that dopants interact only negligibly with 
their periodic images. The bottom layer of the slab was 
passivated; the atomic positions were fully relaxed in all 
but the bottom two layers. 

The wavefunctions and charge density were expanded 
in a plane-wave basis and evaluated using an ultra- 
soft pseudopotential approach2Si2i as implemented in the 
VASP codeiSSiSi^ Electron correlations were treated at 
the level of the local spin-density approximation (LSDA) 
with total energies converged to 10^"* eV. 

A single wave vector was used to sample the Brilloun 
zone of the slab supercells and a plane-wave cutoff of 
227.24 eV was used for all the calculations. For super- 
cells containing Mn, atoms within ~7.0 A of the Mn 
site(s) were fully relaxed until the total energy between 
two structural configurations changed by less than lO"'^ 
eV. 



B. Simulated STM Images 

Theoretical STM images were simulated using the 
method of Tersoff and Hamanni^^ In this approach, the 
central quantity is the the local density of states (LDOS) 
near the surface. The LDOS is integrated over an en- 
ergy range determined by the experimental bias voltage; 
contours of constant energy-integrated LDOS simulate a 
constant-current image. For filled state imaging, the up- 
per bound of this energy range is the Fermi level. The 
Mn-doped GaAs sample under consideration in this work 
is p-type with an estimated hole density of lO^^cm"^, so 
that the Fermi level should be near the valence band 
maximum (VBM). Since we do not know the distribu- 
tion of donors and acceptors present in the sample, this 
information alone is not sufficient to determine the ex- 
act location of the Fermi level for our sample. However, 
we have found that, as a consequence of the large bias, 
the simulated images are rather insensitive to the value 
of the Fermi level used in the LDOS integration, and 
so we choose arbitrarily a point 0.4 eV above the VBM 
as the upper bound for all the STM simulations. For 
this choice, the Fermi level is above the acceptor level 
of substitutional Mn and below the donor levels of inter- 
stitial Mn. Simulations performed with the position of 
the Fermi level shifted by ±0.3 eV give similar results. 
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The location of the VBM for each type of dopant was de- 
termined by inspection of the total DOS relative to that 
projected onto the dopant site. 

As a measure of the agreement between simulated and 
measured STM images, we consider mainly the overall 
shape and spatial extent of features of interest. The 
defect-free (110) surface has a (110) mirror plane sym- 
metry, so it is also useful to note whether a particular 
dopant configuration preserves, even approximately, this 
symmetry. Finally, we use qualitative arguments to judge 
whether the number of observed features are consistent 
with our total-energy calculations. 

IV. SUBSTITUTIONAL MN 
A. Simple Substitutional Dopants 

Substitutional Mn dopants are of primary importance 
for the magnetic and transport properties of Mn-doped 
GaAs, and hence we discuss the simulated STM images 
for these types of dopants first. Since STM is a near- 
surface sensitive probe, we have only considered substi- 
tutional Mn dopants in the top three layers near the (110) 
surface, as shown in Fig.[21[a). 

We will use the notation s{n) to denote a substitutional 
Mn in layer n (with the topmost surface layer defined as 
n = 1). Figures |3Ia), (b), and (c) show the simulated 
images from s(l), s(2), and s(3), respectively. The s(l) 
and s(3) images are strikingly similar, for the following 
reasons. The s(l) Mn atom forms bonds with two As 
atoms that are nominally in the same plane; likewise for 
the s(3) Mn atom. Thus the overall shape of the dopant 
consists of three overlapping maxima, centered on the 
two As atoms and on the Mn atom, giving rise to an 
elongated feature with a (110) mirror-plane symmetry. 
Hence, both of these sites are consistent with feature A 
in Fig.n 

The s(2) dopant on the other hand, has a very dif- 
ferent image. In the s(2) location the Mn forms a bond 
with an As atom in the top layer, with the bond pointing 
toward the surface. The perturbation due to s(2) is cen- 
tered on the surface As site involved in the bond. This 
overlaps with the contribution from the Mn dopant itself, 
and thus the shape is more circular compared to the s(l) 
and s(3) dopants. It is also larger in size than surface As 
atom contributions and has a (110) mirror-plane symme- 
try. These agree qualitatively with the characteristics of 
feature B in Fig. ^ leading us to identify those features 
as arising from s(2) dopants. 

B. Substitutional Mn Pairs 

Since there will always be some fraction of neighbor- 
ing Mn that occur simply by chance, we consider pairs of 
substitutional Mn as well. Figures|2Id), (e) and (f) show 
simulated images of pairs of Mn substitutionals s(l)s(l). 



s(3)s(3) and s(3)s(l), respectively. For the s(3)s(l) pair, 
we note that even though the two Mn involved are '^6.8 
A apart, and in different atomic planes, because of the 
projected view of the STM scan the appearance of this 
pair is quite similar to those of the s(l)s(l) and s(3)s(3) 
pairs. All three of these simulated images have a shape 
comparable to that of feature C in Fig. ^ Common fea- 
tures include a region of low intensity between the two 
dopants, a (110) mirror-plane symmetry (the experimen- 
tal data appear to only approximately have this symme- 
try) and a similar spatial extent involving four surface 
As atoms. It is also interesting to note that these pairs 
have simulated images that are essentially a linear su- 
perposition of the images of the isolated substitutionals 
shown in Figs.Ota) and (c). This suggests that, at the 
level of detail that STM reveals, two substitutional Mn 
do not interact very strongly when they are next-nearest 
neighbors. 



Figures OJg) and (h) show the simulated images 
for nearest-neighbor substitutional pairs s(l)s(2) and 
s(3)s(2), respectively. Overall, the simulated images of 
these pairs are similar to what one would obtain from a 
linear superposition of the images of the isolated defects. 
However, the simulated image for s(l)s(2) clearly reveals 
that the contribution from the s(l) dopant is now asym- 
metric, lacking a (110) mirror plane. This indicates that 
substitutional pairs interact more strongly when they are 
in near-neighbor configuration than in the next-nearest- 
neighbor configuration described in the paragraph above. 
The shape of these simulated images is generally consis- 
tent with that of feature D in Fig. with a character- 
istic "L" shape involving three As atoms, two of which 
are in the same As (001) plane. We have not found any 
other dopants that yield such a favorable comparison to 
the experimental data, leading us to identify feature D as 
s(l)s(2) and s(3)s(2) pairs. In the experimental data one 
can observe three of the four possible crystal orientations 
of these pairs of substitutionals. 



An estimate of the number of neighboring pairs of Mn 
substitutionals based on a random distribution of Mn on 
the Ga sublattice is much lower than the number of com- 
pound features, in particular those labeled "C", seen in 
Fig.n Even if there is some positional correlation among 
the substitutional dopants, for example, the strong short- 
range attraction between Mn found in a recent theoreti- 
cal study,^^ it is unlikely that a description based entirely 
on substitutionals and clusters of substitutionals can de- 
scribe the frequency and shapes of all the features in 
Fig. n This suggests that other Mn complexes, in par- 
ticular those involving interstitial Mn, may help explain 
the remaining features. 
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FIG. 2: (color) Top view (upper) and side view (lower) showing the Mn dopant configurations considered in this work. Ga, As, 
and Mn atoms are shown as cyan, yellow, and red circles, respectively. The dashed lines separate regions containing different 
types of dopants. From left to right: (a) substitutional Mn in the first three layers, (b) interstitial Mn in the first two layers, 
(c) complexes of substitutional and interstitial Mn in layer 2, (d) complexes of substitutional Mn in layer 1 and interstitial Mn 
in layer 2, (e) complexes of substitutional and interstitial Mn both in layer 1. For the complexes in regions (c), (d), and (e), 
several alternative locations of the interstitial Mn involved in the complex are shown in a lighter shade of red. See text for 
explanation of the notation. 



V. DOPANT COMPLEXES 

A. Total Energy Considerations 

It has recently been shown that although the formation 
energy of interstitial Mn is much higher than that of sub- 
stitutional Mn in bulk GaAs, the presence of the (001) 
growth surface significantly enhances the likelihood of 
Mn occupying interstitial sites Thus one expects that 
both substitutional and interstitial Mn will be present 
in typical Mn-doped GaAs samples, as shown by re- 
cent measurements of the Mn distribution in Mn-doped 



GaAsii 

In the bulk, the barrier for diffusion of a charged (+2e) 
interstitial Mn is ~0.5 eV; therefore at room temperature 
we expect interstitial Mn to rapidly diffuse throughout 
the bulk material unless it becomes kinetically trapped. 
Since substitutional dopants are acceptors, they are a 
natural trap for interstitials, electrostatically binding to 
them in a number of different physical configurations. 
Clustering phenomena for compensating dopants such as 
these have recently been predicted using Monte Carlo 
techniques. On the other hand, since interstitial Mn in 
bulk is a double donor with charge of +2e, under condi- 
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FIG. 3: Panels (a), (b), and (c) show show simulated STM 
images of isolated Mn substitutionals. Panels (d) through 
(h) show simulated images for Mn substitutional pairs. Mn 
positions are marked with +. The crystalline orientation is 
the same as Fig. |21 



tions of filled-state imaging (for which the tip is positively 
biased relative to the sample) free interstitial Mn will be 
repelled away from the surface, leaving behind only those 
interstitials that arc bound in complexes. To determine 
which such interstitial-substitutional complexes are most 
likely to be present, we turn to total-energy calculations. 

To simplify the discussion, we extend the notation in- 
troduced earlier to now include Mn interstitials, and com- 
plexes of interstitial and substitutional Mn. Isolated in- 
terstitial Mn in layer n, coordinated to As and Ga, will 
be referred to as i-^^{n) and «*^^(rt), respectively. For ref- 
erence, the configuration of isolated interstitial dopants 
in the first two layers is shown in Fig. El^b) . Complexes 
composed of interstitial and substitutional Mn will be 
referred to with a composite notation that takes into ac- 
count the crystalline orientation of the atoms involved. 
For example, since the GaAs (110) surface does not have 
a (001) mirror-plane symmetry, there are two distinct 
ways of forming a complex of interstitial and substitu- 
tional Mn. In Fig. |21 the surface Ga atoms are to the 
right of the surface As atoms to which they are bonded, 
so that we may consider complexes where the interstitial 
Mn is to the left or right of the neighboring substitutional 
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FIG. 4: Relative total energies of complexes of interstitial and 
substitutional Mn near the GaAs (110) surface: (a) s(l)i°''(l) 
and a{l)i^\l), (b) s(m)i*^"(2) and s(m)i*=(2) (m = 1,2), 
(c) j°''(l)s(l) and i*=(l)s(l), and (d) i°"-(2)s(l). Simulated 
STM images are provided for the low-energy configurations 
in the shaded regions. 



dopant; such complexes will be denoted as i^^in)s{ra) 
and s{m)i^^{n), respectively. Figure EJc), (d), and (e) 
shows the configuration of near-surface complexes we 
have considered in this work. Thus, for example, Fig.jS^c) 
shows the complexes s(2)i^^(2) and s(2)z^''(2), which we 
discuss below. 

Figure^shows the calculated total energies of Mn com- 
plexes relative to that of the s(2)z^^(2) dopant, which we 
identify as the lowest energy configuration. There is a 
large variation in the energies, up to ~1.3 eV above the 
energy of 5(2)1^*^(2). We also note that for complexes in- 
volving only Mn atoms on the surface (Figs.^a) and (c)) 
the lowest energy configuration occurs when the intersti- 
tial site is nearest the substitutional site. Relative to 
these nearest-neighbor complexes, when interstitial Mn 
is displaced by lateral translations in the (110) direction, 
the energy increases dramatically; the next lowest energy 
configurations are ^0.3 eV higher in energy. Thus, diffu- 
sion of surface interstitial Mn from high-energy locations 
(denoted with * and **) to low-energy configurations is 
likely to form nearest-neighbor complexes on the surface. 

It is difficult to estimate, on energetic grounds alone, 
the expected number of such complexes. The presence 
of both the surface and the STM tip create additional 
complications. For example, the energy of interstitial 
Mn depends strongly on its proximity to the GaAs (110) 
surface; an interstitial Mn in the interior of GaAs far 
from the surface can lower its energy as much as '--^l eV 
by moving to a location just below the surface. This 
suggests that interstitial Mn will diffuse from the bulk- 
like region toward the surface; during the STM scans 
these "excess" interstitials will be observed only if they 
are bound in complexes. Such phenomena are difficult to 
quantify and beyond the scope of this work. Instead we 
choose simply to report simulated STM images of several 



6 




FIG. 5: Simulated STM images for Mn substitutional- in- 
terstitial complexes. The locations of substitutional and in- 
terstitial Mn are denote with + and x, respectively. The 
crystalline orientation is the same as in Fig. |21 



low-energy complexes withm the shaded regions of Fig. 01 



B. Simulated STM Images 

1. Substitutional-interstitial complexes 

Figure [S] shows the simulated images for the low- 
energy configurations of substitutional-interstitial com- 
plexes. The simulated image of the lowest energy 
substitutional-interstitial surface complex, s(l)i^^(l), is 
shown in Fig. Ela). There is a reduction in intensity of 
the As surface atom to the right of the s(l) site and in 
the region between the As (001) planes to the left of the 
s(l) site. The simulated images of s(l)i'^''(l)* [Fig.EIc)], 
s(2)iAs(2) [Fig. Eld)], and s{l)i^''{2) [Fig. I^f)] are all 
similar, but the reduction in the simulated intensity be- 
tween the As (001) planes to the left of the s(l) site is 
less apparent. Among these, only the 5(2)1^*^(2) com- 
plex has a (110) mirror-plane symmetry. All four of 
these simulated images are generally consistent with the 
characteristics of feature E in the experimental data of 
Fig. n However, since the s(2)i^^(2) complex is much 
lower in energy ('^0.5 eV) than the other complexes just 
described, and since there appears to be only two such 
features in Fig.Q] we suggest that feature E is most likely 
a s{2)i^^{2) complex. If this is correct, then the ob- 
served apparent displacement of the surface As atom re- 
sults mainly from an electronic effect rather than a real 
atomic displacement. 



The simulated image of the s(l)i^**(l) complex shown 
in Fig. inib) is rather different than those for the com- 
plexes just described. This simulated image obviously 
has a (110) mirror-plane symmetry and a "butterfly" 
shape with a region of lower intensity between the Mn 
atoms. Since there is no (001) mirror plane, the left and 
right wings of the butterfly are clearly different, with the 
contribution from the s(l) Mn more intense than that of 
the Overall, these characteristics are similar to 

those of feature C, which appears to involve four As sur- 
face atoms in neighboring (001) planes. As we discuss be- 
low, there are other energetically competitive complexes 
with similar simulated images, so that we postpone any 
conclusions regarding feature C until all the relevant com- 
plexes have been discussed. 

The simulated image of the s{2)i'^'^{2) complex shown 
in Fig.EJe) does not have a (110) mirror plane symmetry. 
In this configuration there is an oblate region of intensity 
along the As (001) plane and a reduction in the surface 
As atom intensity nearest the interstitial site. On the 
basis of the available STM data, we cannot rule out the 
possibility that some of those features labeled B in Fig.Q] 
result from 5(2)1*^*^(2) complexes. 



2. Interstitial- substitutional complexes 

Figure^lshows the simulated images for the low-energy 
configurations of interstitial-substitutional complexes. 
The simulated images of the i^"(l)s(l), z^"(l)s(l)*, and 
i'^'*(2)s(l) complexes are all very similar in shape, with 
an increase in intensity on the four nearest surface As 
atoms and between the As (001) planes to the left of the 
s(l) site for the first two and to the right of the s(l) site 
for the latter. Obviously, the 1^*^(1)5(1) and i°'^(2)s(l) 
complexes have a (110) mirror-plane symmetry, whereas 
the i^*'(l)s(l)* appears to approximately have this sym- 
metry. The characteristics of these images are similar 
to that of the s(l)«^^(l) complex in Fig.[SJb), but with 
more intensity between the relevant As (001) planes as a 
result of the interstitial location. 

The simulated image of the i*^^(l)s(l) complex shown 
in Fig. IHI^c) has no apparent symmetry. There is an 
increase in intensity between the two surface As atoms 
nearest the «^^(1) site, whereas the s(l) Mn is only very 
weakly visible. Features such as this are not apparent 
in the experimental data, although XSTM scans of lower 
doped samples may help clarify whether they are present. 

From our STM data and the theoretical simula- 
tions, we cannot unambiguously assign feature C to 
one particular complex among the plausible alternatives: 

s(l)iA«(l), i^^{l)s{l), zA«(l)s(l)*, iGa(2)s(l), s(l)s(l), 

s(3)s(3) and s(l)s(3). Considering the large number of 
such features in Fig. ^ and the comparable energies of 
the different complexes, it is likely that more than one 
of these dopants are present in the scan area of Fig. ^ 
Measurements on samples with much smaller Mn concen- 
tration are ongoing and should be able to address these 
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FIG. 6: Simulated STM images for interstitial-substitutional 
complexes. The locations of substitutional and interstitial Mn 
are denote with + and x , respectively. The crystalline orien- 
tation is the same as in Fig.|5| 



issues in more detail. 



VI. SUMMARY 

In summary, we have combined first-principles calcula- 
tions and high-resolution XSTM measurements to char- 
acterize Mn-doped GaAs. XSTM on a (110) cleavage 
plane reveals five reproducible features not found in bulk 
GaAs. Total-energy calculations were used to screen the 
most likely configuration of Mn dopants near the GaAs 
(110) surface. The low-energy configurations were then 
further scrutinized by comparison of simulated and mea- 
sured STM images. These comparisons reveal that there 
are predominantly two types of Mn-related dopants in the 
sample: (1) those involving only Mn substitutional(s), 
including both isolated and pairs of substitutionals, and 
(2) complexes composed of substitutional and interstitial 
Mn. 
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